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We discuss the determination of electromagnetic form factors from the world data on 
electron-nucleus scattering for nuclei Z < 3, with particular emphasis on the derivation of 
the moments required for comparison with measurements from electronic/muonic atoms and 
isotope shifts. 
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I. INTRODUCTION 


The interest in the topic of form factors of light nuclei is threefold: 

1. The wave functions of light nuclei A < 12 can be calculated today with exact methods [jj 
starting from the nucleon-nucleon interaction known from N-N scattering. Electromagnetic form 
factors at large momentum transfer q allow for the most quantitative test of these wave functions 
including their short-range properties. 

2. Modern experiments provide both charge- and matter isotope shifts for many stable and 
unstable nuclei. The radii (rms-radii, Zemach moments) measurable via electron scattering can 
provide the reference to convert shifts to absolute radii. 

3. For the proton there is presently a major discrepancy (~ 0.04/m) between the rms-radii 
determined from electron scattering jij] and electronic Hydrogen Q] versus muonic Hydrogen [4]. 
Similar comparison for light nuclei should provide additional insight. 

In this paper we briefly describe the determination of the most precise electromagnetic form 
factors. Given the recent focus on radii, we will place particular emphasis on the determination 
of the various moments; for the comparison of the form factors with modern theory we refer the 


reader to two review papers 


, 6 ]. 
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II. DETERMINATION OF RADII 


Of particular importance for the determination of the rms-radius R of the nuclear charge density 
p{r ) — which non-relativistically is simply the Fourier transform of the monopole charge form factor 
G c (q) — is the question concerning the range of momentum transfers q that is relevant for the 
extraction of radii. While the answer ’’low q” is the standard one, the question is rarely addressed 
in a quantitative fashion, with the consequence that improper weight is given to particular sets 
of data and new experiments aiming at a determination of R are carried out at uninteresting q' s. 
Using a standard ’’notch-test”, we have studied this question in detail. In fig.l we show the result 
for the particular case of the proton. Only the data 0.5 < q < 1.2 fmT 1 are really sensitive to the 
rms-radius. At lower q, the finite-size effect in the form factor G(q) ~ 1 — q 2 R 2 /6 + ... is too small 
as compared to the uncertainties of the data, at q > 1.2/m” 1 the higher moments such as (r 4 ) 
dominate. Note that for nuclei, which all have rms-radii larger than the one of the proton, the 
maximum sensitivity occurs at correspondingly lower q , i.e. where q 2 R 2 is of similar size. 

The first Zemach moment Q] ( r )( 2 )i which depends on both the magnetic and the charge form 


factor, is needed to predict HFS in atoms. The third Zemach moment 


(r 3 )( 2 ) is needed to relate 


the Lambshift in muonic atoms to the rms-radius. For both moments, the range of sensitivity is 


quite similar to the one for the rms-radius 


EL 



FIG. 1: Sensitivity of the (e,e) data to the lowest moments (r 2 ) and (r 4 ) for the proton. Plotted is the 
change of the moment for a 1% change of the world data in a narrow region A q around a given value of q. 


A second — and rather recent — insight needs to be considered [lOj. The standard procedure 
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of extrapolating G(q) to q = 0, where the rms-radius can be obtained from the slope dG/dq 2 , 
is unreliable. This extrapolation requires a parameterization of G(q) which is fitted to the data. 
Simply taking a convenient form of G(q) can lead to the situation that the corresponding p(r) 
has an unphysical behavior at large radii r, e.g. can be large or not approach zero as we expect 
from our physics understanding of bound systems. Such a behavior at large r can correspond 
to unreasonable moments given the large weight of large r in the moment-determining integrals. 
In g-space this behavior amounts to an unphysical curvature of G(q) which affects the (implicit) 
extrapolation from the (/-region sensitive to R (see above) to q = 0 where R is extracted. 

The only way to avoid this problem is to make sure that the large-r density behaves in a way 
that is consistent with our understanding of the physics of densities. At large r, the density of 
any composite system behaves like a Whittaker function governed by the separation energy of the 
lightest charged constituent 11]. This Whittaker function (with corrections of minor numerical 
impact) can easily be calculated and used to constrain the large-r shape of p(r)\ typically, this 
constraint can be used in the region where p(r) has fallen to less than 1% of its central value. 
Imposing this constraint during a fit of the (e,e) data is practical if parameterizations are chosen 
that have analytical Fourier transforms so that data (in g-space) and constraint (in r-space) can 
be simultaneously fitted. 

The need for constraining the large-r behavior becomes even more apparent once one considers 
the integral / (| R ° p(r) r 4 dr which determines R in the limit i?o —>• oo. In order to get 98% of R 
for e.g. the deuteron, one has to integrate out to 7/m! It is immediately obvious that electron 
scattering cannot in a significant way determine the finite size contribution to the form factor at the 
correspondingly low q of ~ 0.2/m , accordingly, a determination of R to percent-type accuracy 
would be illusionary. The usual ’’solution” for this problem, the use of parameterized form factors, 
introduces an arbitrariness (model dependence) that can be avoided only by introducing the physics 
constraint as discussed above. 


III. DEUTERON 

The determination of form factors for the 1=1 nucleus deuteron is complicated by the need to 
separate the three formfactors Gc, Gm and Gq. Charge (C+Q) and magnetic (M) contributions 
can be separated via forward- and backward-angle data at the same q. Charge monopole (C) and 
quadrupole (Q) contributions can be separated if tensor polarization observables, T 20 in particular, 
are available. For the deuteron the data base is quite extensive, with some 450 data points in the q- 
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region 0.2 < q < 10/m . A fit with the very flexible Sum-Of-Gaussians (SOG) parameterization 


12| of the data (after correction for Coulomb distortion [13]) supplemented by the large-r tail 


constraint yields the form factors with error bars 5G(q) that account for both the random and the 
systematic errors of the data. The detailed comparison of the resulting form factors with modern 
theory has been discussed in [5]. 

The charge rms-radius resulting from this fit is given in table 1 and compared to the preliminary 
value from muonic deuterium, measured by the CREMA collaboration [n|, and the radius that can 


a 


be derived from the accurate knowledge of the triplet n-p scattering length uai . We find excellent 
agreement within the error bar of 0.010/m of the electron scattering result. This agreement is 
particularly relevant with regard to the ” proton radius puzzle” as for the proton the discrepancy 
between (e,e) and pR amounts to a significantly larger ~ 0.04/m. 


(e,e) 

2.130 ±0.010/to 

pR 

2.1289 ± 0.0012/to (prelim.) 

Q"n—p 

2.131/to 


TABLE I: Deuteron charge rms-radii. The pH-value is preliminary due to potential incompleteness of the 
nuclear polarization correction. 


IV. HELIUM 


Electron scattering data are available for 3 He ( 310 data points in the region 0.2 < q < 10/m _1 ) 
and 4 He (190 data points for 0.2 < q < 8.8/m _1 ), the data for 4 He being significantly more precise. 
They are also simpler to analyze as no error-enhancing C/M-separation is needed. 

The Gc and Gm form factors have been determined as described above for the deuteron. A 
detailed comparison to modern theory is discussed in a review paper 5|. Overall, it can be stated 
that the agreement of calculated form factor and experiment is amazingly good, despite the fact 
that the form factors contain substantial contributions due to Meson Exchange Currents MEC. As 
an illustration we show in fig.2 the data for 3 He, together with modern calculations that solve the 
Schrodinger equation for NN potentials derived from NN scattering. 

The rms-radii together with the Zernach moments have been calculated [9!] from these fits and 
are listed in table 2. 

Isotope shifts for the unstable isotopes 6 He and 8 He have also been measured via optical transi- 

nn 

tions in nuclei stored in traps JJJ, |20|], see table 3. For d He several measurements of the 3-4 isotope 
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FIG. 2: Charge and magnetic form factors of 3 He ( poi nts with error bars) together with calculations 
(IA+MEC) of P. Sauer, R. Schiavilla and J. Golak 16 18]. 



3 He 

4 He 

rms-radius charge 

1.973 ± 0.016/m 

1.681 ± 0.004/m 

rms-radius magnetism 

1.976 ± 0.047/m 


{rh 

2.628 ± 0.016/m 


(r 3 h 

28.15 ± 0.70/m 3 

16.73 ± 0.10/m 3 


TABLE II: Moments of Helium isotopes from (e,e). 


shift from optical transitions are available 


21[-[2J], but they scatter by roughly the error bar of 


the electron scattering value. The isotope shift for 3 He is positive due to the more extended proton 
wave function resulting from the lower proton separation energy. Interestingly, the shift for 8 He 
is smaller than the one for 6 He. Due to the more symmetric configuration of the extra neutrons 
- closed p 3 / 2 -shell for 8 He — the center-of-mass movement of the 4 He core in 8 He has a lesser 
amplitude and this leads to a lesser ’’smearing” of the charge. 

Isotope shifts of the matter radii have been deduced via scattering of GeV/nucleon nuclei on 


hydrogen in inverse kinematics 


251 ]. They are also listed in table 3. 



3 (^) charge 

$ Rmatter 

3 He 

1.066 ± 0.060/m 2 


6 He 

1.047 ± 0.034/m 2 

0.81 ± 0.08/m 

8 He 

0.911 ± 0.095/m 2 

0.96 ± 0.08/m 


TABLE III: Charge and matter isotope shifts x He- 4 He relative to 4 He. 


The reference point for all the Helium radii is 4 He. For this nucleus the data from electron 
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scattering is the most precise, see fig.3. In addition, we know for this nucleus not only the shape 
of the large-radius density, but also the absolute value. The extensive set of proton- 4 He elastic 
scattering data has been analyzed using Forward Dispersion Relations FDR llj; the residuum of 


the singularity closest to the physical region, due to exchange scattering at 0°, yields the asymptotic 
norm of the proton wave function to ±10%. This knowledge helps to extract a very precise charge 
radius of 1.681 ± 0.004/m 26] . 



FIG. 3: Ratio of data/fit for 4 He at low q. The data are from refs. 



This 4 He charge rms-radius actually is the most precise radius for any nucleus determined by 
electron scattering. It therefore is particularly interesting to compare this radius to the one being 


determined from muonic Helium 


311 ]. The final value from this experiment is not yet available, 


but it is known to be well within the error bar of ±0.004/m of the radius from electron scattering. 
This agreement of better than 0.004/m is highly significant given the 0.04/m discrepancy between 
(e,e) and pX for the proton. This speaks strongly against a potential difference between the 
’’electromagnetic” interaction of electrons and muons, an idea that had been speculated about in 
order to explain the proton radius puzzle. 


V. LITHIUM 

For Lithium electron scattering data are available for A=6,7, the latter being less accurate and 
less extensive; moreover they contain both contributions from the unresolved scattering to the 
first excited state of "Li and a large contribution from quadrupole scattering. 6 Li is therefore the 
natural reference nucleus. 
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321 ]; the 


The 6 Li charge form factor and rms-radius has been determined as discussed above 
data set comprises 86 points in the range 0.1 < q < 3.8/m” 1 . A potential complication arises 
due to the partial a + d cluster structure and the low deuteron removal energy. For this case, the 
asymptotic behavior may be more complex. For the shape of the tail we have therefore used the one 


from the GFMC-calculation of Pieper et a/.jll!]. As this calculation corresponds to an exact solution 
of the Schrodinger equation and reproduces the experimental binding energy, the large-r shape can 
be trusted. The resulting charge rms-radius amounts to 2.589 ± 0.039/m. The comparatively large 
uncertainty results from the limitations in the (e,e) data. 


33] while for the Ml 


The charge form factor for 6 Li agrees very well with the GFMC calculation 
form factor, which also has been measured up to q ~ 2.8/m _1 , a small disagreement is observed 
at the highest g’s. 

Isotope shifts of the charge rms-radii have been measured by Nortershauser et al. (for references 


see 


321 ]) for mass number up to A=ll. n Li is a particularly interesting case of a Borromean nucleus, 


as all sub-systems, 9 Li+n and n+n (di-neutron), are unbound while the 3-body system 9 Li+n+n 
is bound. These shifts have been measured at CERN using laser spectroscopy with the Li-nuclei 
stored in traps. Fig.4 shows the resulting radii, together with a series of theoretical calculations 
discussed in detail in [32|]. The unusual behavior of the radii is related to the fact that 6 Li and ' Li 


show pronounced effects of cluster structure while 8 Li and 9 Li are closer to a mean-field description. 



FIG. 4: Charge radii of Li-isotopes (reference point 6 Li) from 


321 ]. The data are from refs. 
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The matter radii of the Li-isotopes have been measured at GSI [37(. Experimentally these 
quantities are accessible (although with lesser precision) by scattering of Lithium nuclei, produced 
by fragmentation with energies of the order GeU/nucleon, from Hydrogen; cross sections in this 
energy range can be interpreted using Glauber theory. The resulting data show a more or less 
constant matter radius throughout the series A=6-9, with a 30% increase for A=ll resulting from 
the extremely low two-neutron separation energy of only 370 keV. 


VI. MAGNETIC FORM FACTORS 


Magnetic scattering for most nuclei is dominated by the properties of valence neutrons or 
protons, thus providing an observable governed by different physics. In general, magnetic form 
factors also receive larger contributions from meson exchange currents which, depending on the 
research goal, is a benefit or a complication. 

We have in the discussion given above placed little emphasis on these magnetic form factors. 
Magnetic form factors are more difficult to determine as they require experiments at very large 


scattering angle. As a consequence, t 


form factors has been published in 


re data base is more limited. A review paper on magnetic 


38|. 


Magnetic rms-radii are particularly difficult to determine as at low q the contributions from 
charge scattering dominate the electron scattering cross sections. Charge scattering can be sup¬ 
pressed by measurements at 180° scattering angle; since a number of years, however, no such 180° 
facilities are available anymore. 


VII. SUMMARY 

Precise rms-radii from electron scattering are of interest for comparison with radii obtained from 
modern ab initio calculations and measurements from muonic X-ray data; they also serve as an 
anchor point for the many isotope shifts now measurable for unstable nuclei. In this paper, we have 
used the world data on electron scattering to determine accurate radii and Zemach moments for 
Z < 3. The uncertainties of the radii derived from (e,e) are between 0.5-1.5%. For the special case 
of 4 He, where the highest accuracy of ~0.25% is reached, we find good agreement with the radius 
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from muonic Helium, contrary to the case of the proton 


39] where a 4% discrepancy remains. 
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